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1.0  INTRODUCTION 

In  this  report,  three  widely  used  methods  for  imposing  discrete  solid-wall  boundary 
conditions  for  the  numerical  computation  of  inviscid  flows  are  presented,  and  two  of  these 
methods  are  shown  to  offer  only  provisionally  correct  results.  Specifically,  it  is  shown  that 
a  method  involving  contravariant  velocity  components  can  generate  distortions  in  computed 
flows.  This  method  has  been  used  for  years  in  the  AIR3D  code  (Ref.  1)  and  its  derivatives 
to  compute  steady-state  flows  in  stationary  grid  systems.  Also,  with  the  inclusion  of  time 
metrics,  this  method  has  been  used  to  compute  flows  in  moving  grid  systems  (Ref.  2). 

Steady-state  computations  performed  with  these  codes  have  shown  that  distortions  can 
be  generated  in  computed  flows  when  grid  lines  emanating  from  a  solid  boundary  are  not 
orthogonal  to  the  wall.  The  consequences  of  grid  skewness  are  discussed  in  Ref.  3.  In  addition, 
it  is  shown  here  that  the  method  of  contravariant  velocity  components  is  correct  only  if  grid 
lines  emanating  from  solid  walls  are  normally  oriented.  Because  complex  CFD  calculations 
can  be  difficult  to  interpret,  the  claims  made  here  are  demonstrated  with  explicit  solutions 
to  simple  model  problems. 

Instead  of  the  above  procedure  for  determining  velocity  at  solid  boundaries,  there  is  a 
more  direct  approach  of  simply  extrapolating  velocity  to  the  wall  and  subtracting  the  normal 
component.  Such  a  method  gives  correct  results  on  a  stationary  grid  without  conditions  on 
the  grid,  such  as  the  orthogonality  constraint  mentioned  above.  While  this  more  direct 
approach  has  an  unambiguous  implementation  for  stationary  grid  systems,  there  are  various 
ways  to  employ  the  idea  for  moving  grid  systems.  Furthermore,  some  of  these  implementations 
for  dynamic  grid  systems  can  cause  spurious  features  in  computed  flows.  A  correct 
implementation  is  proposed  and  compared  v/ith  another  apparently  natural  one.  Again,  since 
complex  CFD  calculations  can  be  difficult  to  interpret,  comparisons  are  made  on  the  basis 
of  explicit  solutions  to  simple  model  problems.  Finally,  the  correct  method  is  shown  to  be 
consistent  with  solid-wall  boundary  conditions  imposed  in  Ref.  4.  Also,  an  explanation 
is  given  for  the  use  in  Ref.  5  of  a  solid-wall  boundary  condition  for  pressure  in  moving 
grid  systems. 

The  contents  of  the  report  can  be  summarized  as  follows.  In  Section  2.0,  the  necessary 
notations  and  assumptions  are  presented.  In  Section  3.0,  definitions  are  given  for  three  methods 
of  imposing  discrete  solid-wall  boundary  conditions  in  moving  grid  systems  for  the  numerical 
computation  of  inviscid  flows.  As  appropriate,  the  boundary  conditions  given  are  simplified 
for  the  treatment  of  stationary  grid  systems.  Then  in  Section  4.0,  simple  examples  are  presented 
to  demonstrate  which  boundary  conditions  from  Section  3.0  give  correct  numerical  results. 
In  Section  5.0,  the  numerical  boundary  conditions  used  in  Refs.  4  and  5  are  discussed.  Finally, 
in  Section  6.0,  the  conclusions  of  this  study  are  summarized. 
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2.0  NOTATIONS  AND  ASSUMPTIONS 

Before  any  discussion  of  discrete  boundary  conditions  can  begin,  the  continuum  boundary 
condition  must  be  given  together  with  a  description  of  the  kind  of  applications  for  which 
it  is  intended.  Also,  since  CFD  codes  do  not  all  operate  in  the  same  way,  the  present  numerical 
orientation  must  be  explained  and  certain  notations  must  be  developed. 

Assume  that  a  solid  body  is  in  rigid  motion  with  respect  to  what  will  be  called  a  stationary 
frame.  Let  stationary  frame  Cartesian  coordinates  be  given  by 

*  =  [x,y,zf  =  [xi,X2,X3f . 


Also,  let  curvilinear  coordinates 

€  =  =  [6,6, 

be  established  around  the  moving  body  for  the  purpose  of  defining  a  computational  grid. 
Specifically,  at  every  time  level,  grid  points  f)  satisfy 

j  k  =  i  =  C{^'jkhn-  (1) 


That  is,  the  curvilinear  coordinates  assume  integer  values  at  grid  points.  Also,  it  is  assumed 
that  for  every  fixed  time,  t,  the  transformation  from  ®  to  ^  is  well-behaved  in  the  sense  that 


6  6  6 

Vx  Vy 

Cx  Cy  Cz 


satisfies 


0  <  7i  <  det  {J}  <  72  <  oo 


(2) 

(3) 


for  constants  7j  and  Therefore,  the  inverse  transformation  from  f  to  x  is  well-behaved. 
Also,  it  follows  from  the  chain  rule  expressions. 


dxj  d^k  _  'V~'  ^6  ^^k  _  c 

^  d^k  dxj  ~  dxk  d^j  “ 


(4) 
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that  the  inverse  of  J  is  given  by 


Xg 

H 

Vi 

Vr, 

y< 

(5) 


Let  the  body  surface  be  defined  by 

C{x,t)  =  constant, 


and  suppose  the  points  away  from  the  body  are  associated  with  increasing  f .  Therefore,  for 
fixed  time,  t,  an  outwardly  directed  unit  normal  vector  at  the  body  surface  is  given  by 


n{x,t) 


VC 

ii'^cir 


(6) 


Next,  let  represent  the  Cartesian  components,  relative  to  the  stationary  frame,  of 
the  trajectory  of  a  material  particle.  It  is  assumed  that  the  curvilinear  coordinates  associated 
with  remain  constant.  Thus,  the  grid  coordinate  system  does  not  move  in  relation  to 
the  body  as  would  be  the  case,  for  example,  if  the  grid  were  adapting  to  the  flow.  In  addition, 
let  x^(t)  represent  the  Cartesian  components,  relative  to  the  stationary  frame,  of  the  trajectory 
of  a  grid  point,  with  necessarily  fixed  curvilinear  coordinates.  Suppose  that  the  curvilinear 
coordinates  associated  with  and  are  the  same, 

i{xm(t),t)  =  ^{xg(t),t)  =  constant.  (7) 


uriiiWj  Uy  LiiC  Li  litcLllv/iX  liwiil 


^  t 

tv/  ^  ■ 


tr  CI>\T£>V\T 

lO  llJl  V  VI  tit/AW  AV/l  vVvijr  |,liAAV)  Aj  IL  AV/lli-7vVo 


that  x(t) 


Thus,  ®„,(0 


where  v  =  dx  /dt  and  v  =  dx Jdt  represent. 


-gv-/,  mm  g  g 

respectively,  the  Cartesian  components,  relative  to  the  stationary  frame,  of  the  velocities  of 
a  material  particle  and  of  a  grid  point. 


To  ensure  that  the  grid  is  in  rigid  motion,  and  thus  in  the  same  rigid  motion  as  the  body 
according  to  the  above  discussion,  the  following  assumption  is  made.  Let  two  points  be  chosen 
arbitrarily  with  fixed  curvilinear  coordinates  and  let  xj^t)  and  *^(t)  denote  the  Cartesian 
components  of  their  trajectories  relative  to  the  stationary  frame.  Then  it  is  required  that  the 
distance  between  these  two  points  remain  constant,  i.e., 

W^bit)  -  *a(<)||  =  11*6(0)  -  *a(0)||. 
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As  shown  in  Appendix  A,  this  condition  has  the  important  consequence  that  there  is  an 
orthogonal  rotation  matrix,  L{t),  i.e., 

L-^{t)  =  L^{t)  (8) 

such  that 

[xk{t)  -  =  Lit)  [*fr(0)  -  *,(0)] .  (9) 

In  particular,  let  one  point  correspond  to  any  grid  point  and  let  the  other  correspond  to  a 
reference  point  such  as  the  body’s  center  of  gravity.  Thus,  grid  movement  can  be  accomplished 
knowing  only  the  rotation  matrix  and  the  trajectory  of  a  reference  point. 

The  physical  condition  to  be  imposed  on  velocity  for  a  solid-wall  boundary  condition 
in  inviscid  flow  is  that  no  fluid  pass  into  the  body.  Mathematically,  this  means  that  the 
components  of  fluid  velocity  relative  to  the  moving  frame  can  have  no  normal  component 
at  the  body  surface.  It  is  shown  in  Appendix  B  that  this  condition  can  be  stated  in  the  stationary 
frame  as  follows, 

n-  (v  -  =  0  (10) 

where  v  =  [m,v,m']^  represents  the  Cartesian  components  of  fluid  velocity  relative  to  the 
stationary  frame.  For  a  more  computationally  convenient  version  of  the  above  equation, 
the  material  particle  velocity,  can  be  replaced  with  the  grid  velocity,  v^,  as  explained 
above.  Thus,  the  continuum  condition  for  fluid  velocity  can  be  expressed  as 

n-(®-»J  =  0.  (11) 


Also,  can  be  calculated  explicitly  by  differentiating  Eq.  (7)  with  respect  to  t  to  obtain 

o 


^  +  ^  =  0 
dt  ^  dt 


Recalling  that  J  -  the  last  equation  can  be  written  as.7®  +  =  0,  or  by  Eq.  (3), 
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Now,  combining  Eqs.  (6),  (12),  and  (11)  gives  the  condition 


VC  „  ,  VC 

livcil  livcil 


or 


u 

r 

Xr^  X(^ 

'6  ■ 

V 

+ 

a  Cy  G 

Vi  Vn  yc. 

nt 

w 

.  2:^  Zr^  Z(^  _ 

. 

Finally,  using  Eq.  (4),  the  continuum  condition  can  be  written  as 


UCcO  +  VCy  +  W(z  +0  =  0. 

3.0  DEFINITION  OF  METHODS 


(13) 


Equation  (11)  above  shows  that,  for  a  grid  point  on  the  body  surface,  the  normal 
component  of  fluid  velocity  must  be  set  equal  to  that  of  grid  velocity.  However,  this  condition 
gives  no  indication  of  how  the  tangential  components  should  be  calculated.  This  is  the  crucial 
point;  after  all,  there  are  many  ways  to  satisfy  Eq.  (11),  including  setting  «  =  To 
guarantee  continuity  of  the  tangential  component  at  the  wall,  some  form  of  extrapolation 
from  the  field  must  be  used.  First,  a  method  involving  contravariant  velocity  components 
is  described.  Then,  two  other  more  direct  methods  are  defined. 

3.1  METHOD  1 


For  a  description  of  Method  1,  a  careful  definition  is  given  first  for  the  contravariant 
components  of  vectors  in  four-dimensional  space-time.  Let  ®^(t)  =  [Xj{t),  yj(t),  Zj(t)]^ 
represent  the  Cartesian  components,  relative  to  the  stationary  frame,  of  the  trajectory  of 
a  fluid  particle,  so  that  its  world  line  is  represented  by  [x^(t),  Zf{t),  t\^.  In  other  words. 


with  hy,  e^,  denoting  a  Cartesian  basis  for  space-time,  the  world  line  is  given  by 


IV  =  +  yfCy  +  zjk:,  +  tet- 


Its  temporal  derivative  is  given  by 

F  =  uCx  +  vcy  +  wig  +  ef 
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where  again  v  =  [UjV.vv]^  represents  the  Cartesian  components  of  fluid  velocity  relative  to 
the  stationary  frame.  Now  suppose  curvilinear  coordinates  (l.tj.f.r)  = 
established  for  space-time.  Also,  suppose  that  Xi^,ri,^,T)  is  defined  so  that  curvilinear 
coordinate  curves  are  obtained  by  fixing  three  of  the  arguments  in  X  while  varying  the  fourth. 
Thus,  a  principal  basis  for  space-time  is  given  by  {X^,  X^,  X^,  X^},  and  the  corresponding 
reciprocal  basis  is  given  by  {2>^,  Pf ,  Vt)  where,  for  example,  is  defined  by 


"I"  "I"  "I" 


(14) 


These  bases  have  the  property  that^^.  •  V^j  =  6^.  Also,  X^.  is  tangent  to  the  coordinate 
line  with  constant  while  V^j  is  orthogonal  to  the  hypersurface  of  constant  For  a 
4-vector,  say  U,  the  covariant  components,  (Wj,Z/2»^3’^4)’  defined  by 

4  4 

U  =  where  li  ■X^.  =  =  Hi. 

j=i 

Also,  the  contravariant  components,  (ll^,U^,W,U‘^),  are  defined  by* 

4  4 

j=i  i=i 


Therefore,  the  contravariant  components,  {U,V,W,S},  of  V  satisfy 

V  =  UX^  +  VXr,  +  WX^  +  SXr 


where 

U  =  V-V^ 
V  =  V-Vt) 
W  =  V  VC 
S  =  V  Vt 


uCx  +  vCy  -I-  wCz  +  Ct 

uri^  +  vr)y  -I-  WTfz  +  Tft 

uCx  +  vCy  +  wCz  +  C< 

UTx  +  VTy  +  WTz  -f  Tj.  (15) 


*  This  terminology  originates  from  the  fact  that  if  the  coordinate  system  were  changed,  the 
transformation  from  the  old  to  the  new  covariant  components  has  the  same  form  as  the  transformation 
from  the  old  to  the  new  principal  basis  vectors.  On  the  other  hand,  the  transformation  of  the 
contravariant  components  is  different  from  that  of  the  principal  basis  vectors,  but  has  the  same  form 
as  that  for  reciprocal  basis  vectors. 
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Finally,  in  view  of  Eq.  (13),  the  continuum  boundary  condition  can  be  expressed  as 


W  =  0 


(16) 


where  W  is  the  contravariant  component  of  the  space-time  vector  V  given  above. 

It  has  become  customary  to  refer  to  U,  V,  and  W  as  the  contravariant  components  of 
velocity  even  when  $,  rj,  and  f  are  time-varying.  However,  this  use  of  terminology  from  tensor 
analysis  is  not  strictly  correct.  In  fact,  by  Eq.  (3),  {x^,  and  {V|,Vij,Vf}  form  principal 

and  reciprocal  bases  for  three-dimensional  space  at  any  given  time  t.  Calculations  similar 
to  those  prior  to  Eq.  (15)  show  that  the  contravariant  components  of  velocity  are  w  •  V^, 
V  •  Vt?,  and  V  •  Vf,  i.e.. 


»  =  (w  •  +  (®  •  V77)x^  -I-  (»  •  VC)x^  =  {u  -  6)*^  +  {v  -  rit)xj,  +  {W  -  Ct)x^.  (17) 

Of  course,  these  components  coincide  with  U,  V,  and  W  when  r},  and  f  are  time-independent. 

The  solid- wall  boundary  condition  as  shown  in  Eq.  (16)  has  been  implemented  numerically 
according  to  a  procedure  that  is  distilled  below  to  a  single  equation,  Eq.  (21).  This  equation 
is  not  used  widely,  but  it  is  important  for  comparison  with  other  methods.  Also,  it  provides 
a  means  for  representing  the  present  method  in  compact  form.  There  are  many  variants  of 
this  method,  but  the  condition  derived  serves  as  a  representative  of  procedures  involving 
contravariant  components. 

The  numerical  implementation  of  Eq.  (16)  proceeds  as  follows.  First,  the  Cartesian 
components  of  fluid  velocity  are  extrapolated  to  the  solid  surface,  where  C  =  /  =  1 ,  from  the 
first  point  off  the  wall,  where  f  =  /  =  2.  Then  the  components,  U,  V,  and  W are  determined. 
By  Eqs.  (15)  and  (2)  this  process  can  be  written  as 


■  u  ' 

n 

6 

n 

U 

y 

zz 

Vt 

+ 

V 

w 

Ct 

w 

jkl 

jkl 

(18) 


t  If  T  is  independent  of  x,  y,  and  z,  then  the  vector  sets,  and  {X^,X^,X^},  are  identical. 
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Following  Eq.  (16),  the  next  step  is  to  set  fV  to  zero  by  multiplying  the  last  equation  by 
[I  -Ml,  where 


M  = 


0  0  0 
0  0  0 
0  0  1 


(19) 


Next,  the  result  of  this  multiplication  is  translated  to  Cartesian  components.  Again,  using 
Eqs.  (15)  and  (2),  this  process  can  be  written  as. 


u 

V 

w 


n 


jkl 


6 

m 

Ct 


n 

+  [I-M] 

jkl 


u 

V 

w 


n 


jkl 


(20) 


Now,  ^  is  explicitly  equated  withw^^j  and  Eq.  (18)  is  substituted  into  Eq.  (20)  to  obtain. 


6  ' 

n 

/  \  —I 

\ 

6 

n 

U 

n 

m 

1 

1 

+  ^jkl 

V 

> 

Ct 

jkl 

i 

.  Ct  . 

jkl 

W 

jk2> 

As  seen  below,  the  expression  in  braces  can  be  expressed  in  terms  of  fluid  velocity  and  grid 
velocity  by  factoring  out  For  this,  the  product,  inserted  after  [/  -M']. 

After  some  algebra,  the  following  is  obtained, 


The  bracketed  expression  can  be  simplified  with  the  following  explicit  calculation. 


■  ■ 

n 

■  0 

0 

0  ■ 

I- 

•  •  VC 

0 

0 

0 

jkl 

0 

0 

1 

1 - 

N 

1 _ 

1  ■ 
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Combining  the  last  two  equations  and  Eq.  (12)  gives 


-n+l  _ 


=  ('»)?«  + 


+  ^Jta} 


or 


n 

jk2 


Mm  ■ 


(21) 


This  will  be  referred  to  as  Method  1  for  setting  inviscid  flow  velocity  at  a  solid-wall  boundary. 
The  crucial  element  on  which  to  focus  is  the  presence  of  the  term,  This  vector  is  tangent 
to  the  coordinate  line  emanating  from  the  the  body  surface,  and  it  is  not  normally  oriented 
unless  and  Vf  are  parallel.  In  Subsection  4.3,  it  is  shown  that  this  aspect  of  Method  1 
causes  spurious  disturbances  to  be  generated  in  computed  flows. 


Finally,  to  complete  the  definition  of  Method  1,  the  grid  velocity  is  specified  with  a  discrete 
formula  instead  of  Eq.  (12).  Since  t>^  =  dx^/dt,  the  local  grid  velocity  at  time  level,  f,  can 
be  approximated  using  grid  point  coordinates  at  time  levels,  f  and 


(„  Y  -  3kl 

y^9)iki  - 


At 


(22) 


Now,  instead  of  storing  the  old  grid,®^^y*  can  be  computed  from  as  suggested  following 
Eq.  (9).  Specifically,  suppose  the  grid  movement  is  characterized  by  the  trajectory  of  a  reference 
point  XQ(t)  and  a  rotation  matrix,  L(t).  Then  by  Eq.  (9), 


Hn 


and 


^°3kl 


Eo(t°) 


Using  Eq.  (8)  to  eliminate  —  *Q(t°)]  between  these  equations  gives 


Xjki  ®  0 


(23) 
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Thus,  by  combining  Eqs.  (22)  and  (23),  the  grid  velocity  is  computed  equivalently  according 
to  Eq.  (22)  or 


(”.)"«  =  i;  {[*”«  -  *««“■')]  - 


-  »o(i”)] }  ■ 


3.2  METHOD  2 


(24) 


Method  2  for  setting  the  inviscid  flow  velocity  is  a  more  direct  numerical  formulation 
of  Eq.  (1 1),  which  has  also  been  used  in  AIR3D.  The  basis  of  the  method  is  as  follows.  Working 
in  the  moving  frame,  the  velocity  at  time  level,  f,  is  extrapolated  to  the  solid  surface  from 
the  first  point  off  the  wall.  Then  the  normal  component  is  subtracted,  and  the  result  is  expressed 
in  the  stationary  frame.  Finally,  the  velocity  obtained  is  used  at  the  wall  for  time  level, 
Using  the  techniques  of  Appendix  B,  and  Eqs.  (B-6)  and  (B-7)  in  particular,  this  procedure 
can  be  expressed  purely  in  the  stationary  frame  as 


•IS  =  (•»)"«  + 


jkl 


I  —  hn^  (w  —  Va) 


n 

jk2  ' 


Using  Eq.  (6),  Method  2  can  be  written  as 


+ (•  - -  jj^p  ■  ('  -  •»)”« 


(25) 


with  the  grid  velocity  computed  according  to  Eq.  (24). 

Method  2  appears  natural  because  it  is  consistent  with  a  standard  condition  posed  in  the 
moving  frame  with  respect  to  which  the  grid  velocity  is  zero.  However,  as  shown  in  Subsection 
4.4,  Method  2  causes  spurious  disturbances  to  emerge  in  computed  flows  because  is 
evaluated  at  more  than  one  point  in  the  above  equation.  Specifically,  this  creates  a  problem 
in  cases  where  the  moving  frame  is  noninertial  and  the  grid  velocity  varies  throughout  the  grid. 

3.3  METHOD  3 


Finally,  Method  3  for  setting  the  inviscid  flow  velocity  is  derived  from  Method  2  by  making 
modifications  to  avoid  the  spurious  effects  shown  in  Subsection  4.4  in  connection  with  Method 
2.  Specifically,  the  present  method  is  posed  so  that  the  grid  velocity  is  always  evaluated  at 
a  single  point.  Method  3  can  be  written  compactly  as 


_n+l  _ 
^jkl 


'jkl 
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Using  Eq.  (6),  this  can  be  written  as 


71+1  _  _ ]  \ 

^jk\  ~  ^jk2  [^Jk2 


(26) 


with  the  grid  velocity  computed  according  to  Eq.  (24).  It  is  shown  in  Subsection  4.5  that 
Method  3  passes  the  tests  developed  in  the  next  section. 


4.0  ANALYSIS  OF  METHODS 

Methods  1 ,  2,  and  3  are  analyzed  in  terms  of  their  performance  for  two  simple  model 
problems.  In  the  first  model  problem,  the  inviscid  flow  is  over  a  flat  plate.  Also,  the  grid 
is  stationary  and  grid  lines  emanating  from  the  plate  have  a  nonnormal  orientation.  In  the 
second  example,  the  inviscid  fluid  is  at  rest  around  a  rotating  cylinder.  Here  the  grid  is  in 
rigid  motion  with  the  solid  body. 


4.1  PROBLEM  1 


For  model  Problem  1,  let  a  flat  plate  be  located  at  z  =  0  and  the  field  at  z  >  0.  Also, 
assume  that  the  fluid  velocity  distribution,  ®  =  [u,v,w\^,  is  two-dimensional  with  v  =  0.  The 
components,  u  and  w,  can  be  chosen  arbitrarily  to  satisfy  the  Euler  equations,  provided 
w  =  0  only  at  z  =  0.  Thus,  a  proper  boundary  condition  will  not  introduce  a  nonzero  v 
component  to  the  velocity. 

The  curvilinear  coordinate  system  for  Problem  1  is  time-independent  and  is  defined  as 
follows, 

^  =  a;,  ■q  =  y-'^z,  (  =  z.  (27) 

To  show  the  skewness  of  this  coordinate  system,  the  orientation  of  the  coordinate  axes  is 
now  computed.  For  example,  the  ^-axis  is  defined  as  the  line  along  which  rj  and  f  are  constant. 
So,  the  position  vector,  *  =  {x,y,z]^,  is  a  function  only  of  ^  along  this  line.  Therefore, 
is  tangent  to  the  ^-axis.  Similarly,  and  give  the  orientations  of  the  ij-  and  f-axes, 
respectively.  These  can  be  calculated  from 

*  =  y^v  +  iC,  ^  =  C- 
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Specifically,  the  orientations  of  the  curvilinear  coordinate  system  axes  are  given  by  the  principal 
basis  vectors  (see  Fig.  1), 


—  “y  ”1”  • 


Note  that  the  and  jj-axes  are  aligned  with  the  x-  and  j-axes,  respectively.  However,  the 
f-axis  lies  in  the^,z  plane.  Not  only  is  it  not  aligned  with  the  z-axis,  it  approaches  the^axis 
as  7  -*  00.  Finally,  note  that  the  reciprocal  basis  vectors  are  given  by  (see  Fig.  2), 

=  Cx,  Vrj=€y-'yez,  VC  =  e^. 


Figure  2.  Reciprocal  basis  vectors  for  Problem  1. 

These  vectors  are  orthogonal  to  surfaces  of  constant  C,  t),  and  f ,  respectively.  Also,  with 
Eq.  (4),  recall  that  a  given  reciprocal  basis  vector,  is  orthogonal  to  the  two  principal 
basis  vectors,  and  has  unit  scalar  product  with  its  counterpart,  x^,. 

4.2  PROBLEM  2 

For  model  Problem  2,  let  a  rotating  cylinder  be  located  at  :>^  +  -  1,  and  the  field 

at  >  1.  Also,  assume  that  the  fluid  velocity  distribution  satisfies  »  =  [u,v,w]^ 

=  0.  Thus,  a  proper  boundary  condition  will  not  introduce  nonzero  components  to  the  velocity. 
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The  curvilinear  coordinate  system  for  Problem  2  is  time- varying  and  is  defined  as  follows, 


e  =  {N-l)z  +  l 

T]  =  1  +  (iV  -  1)  t  -  tan"^(j//x)] 

C  =  {N  -  +  2/2  _  (jv  _  2) 


T  =  t. 


(28) 


Here,  Nis  a  grid  refinement  parameter,  i.e.,  the  approximate  grid  cell  volume  is  1/(A'— 1)^. 
Recall  that  the  cylinder  is  rotating  in  rigid  motion  with  this  coordinate  system,  i.e.,  any  material 
particle  holds  fixed  curvilinear  coordinates.  The  Cartesian  coordinates  can  be  written  in  terms 
of  the  curvilinear  coordinates  as  follows. 


X 


<  +  {N-2) 
(N-1) 


COS 


(^- 1)' 
(iV-l)J 


y 


<  +  (N-2) 
(N-l) 


sin 


(^-1)1 

(iv-i)J 


(^-1) 

(N-l) 


t  =  r. 


(29) 


Thus,  the  space-time  principal  basis  vectors  are  (see  Fig.  3) 

„  — 


C  +  (iV-2)  . 

ClTl 

\t  1 

_ 

C  +  (Ar-2) 

L  (»?- 

1)1 

{N-lf 

(iV-l)J 

{N-lf 

{N- 

■1)J 

(N-1)^  (N-l) 


1 


(N-l) 


COS 


T  — 


in-i) 


(N-n 


ex  + 


(N-l) 


sm 


T  — 


(»?- 1) 


{N  -  1) 


-tK  + 


{N  -  l)^^M^  ^  '  (iV  -  l)^^2TF 
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C  +  (iV-2) 
{N-l) 


sin 


in  - 1) ' 

{N-l)\ 


^  I 

ex  + 


C  +  {N-2) 
(N-l) 


COS 


(^- 1) ' 
(N  -  1)J 


ey  +  et 


-ye^  +  X€y  +  et 


Figure  3.  Principal  basis  vectors  for  Problem  2  (z  and  ^  coordinates  suppressed). 

and  the  reciprocal  basis  vectors  are  (see  Fig.  4) 

=  {N-l)e, 


Figure  4.  Reciprocal  basis  vectors  for  Problem  2  (z  and  ^  coordinates  suppressed). 
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4.3  ANALYSIS  OF  METHOD  1 

Now,  consider  the  application  of  Method  1  in  Eq.  (21)  to  Problem  1 .  Since  the  curvilinear 

coordinate  system  is  time-independent,  the  grid  velocity  is  zero.  Also,  since  x.  =  7e  +  e 

i  y  "z 

and  Vf  =  e^,  Eq.  (21)  becomes 


for  the  first  time  step.  Note  that  the  last  equality  follows  since  it  is  assumed  that  the  flow 
is  two-dimensional  with  the  v  component  zero.  However,  this  method  creates  a  v  component, 
vj^j  =  -  71^9^2’  the  wall  which  is  propagated  into  the  field.  In  fact,  the  problem  is  made 
worse  with  more  skewness  in  the  grid,  i.e.,  as  T  -*  «>.  Since  Method  1  generates  this  disturbance 
in  the  flow,  it  is  judged  to  be  unacceptable. 

On  the  other  hand,  if  the  nontangential  grid  lines  are  orthogonal  to  the  body  surface. 
Method  1  produces  correct  boundary  values.  This  can  be  seen  by  setting  y  =  0  in  Problem 
1  so  that  and  vj^j  =  0.  In  general,  under  the  condition  that  nontangential  grid  lines 

are  orthogonal  to  the  body  surface.  Method  1  can  be  rewritten  as  follows.  Equation  (21) 
must  be  modified  to  reflect  the  constraint  that  must  be  orthogonal  to  the  wall.  Since  Vf 
is  always  orthogonal  to  the  wall,  it  must  be  that  =  cYf  for  some  scalar  c.  On  the  other 
hand,  since  and  Vf  are  principal  and  reciprocal  basis  vectors,  respectively,  from  Eq.  (4) 
it  follows  that  x^.  •  Vf  =  1.  Therefore, 

1  =  .  VC  =  cvc  ■  VC  =»  *<  =  • 

Under  the  new  constraint.  Method  1  takes  the  form, 

-S'  =  »”«  -  jj^p  {VCS:.  •  [.^2  -  (.,&.]  } 

with  Wg  computed  according  to  Eq.  (24).  In  other  words,  under  the  orthogonality  condition, 
Eqs.  (21)  and  (26)  are  identical. 

The  flaw  in  Method  1  can  be  understood  most  easily  as  follows.  For  simplicity,  assume 
that  the  grid  system  is  time-invariant.  Then  the  fluid  velocity  can  be  written  as  shown  in 
Eq.  (17), 

B  =  (®  •  VC)«^  +  (w  •  V?/)x,j  -f  (w  •  VC)x<  =  17*4  +  +  IT *4. 
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Consider  this  in  relation  to  Fig,  1.  Since  and  always  lie  in  the  tangent  plane  of  the  body 
surface,  setting  ®  •  Vf  =  W  =  0  does  indeed  make  v  purely  tangential.  However,  with  skewness 
in  the  grid,  has  a  tangential  part.  Thus,  setting  W  =  0  modifies  the  tangential  part  of  v. 

4.4  ANALYSIS  OF  METHOD  2 

Now  consider  the  application  of  Method  2  in  Eq.  (25)  to  Problem  2.  Since  the  fluid  velocity 
is  assumed  to  be  zero  everywhere  initially,  the  method  becomes 


\^9)jk2  |j  |P 


(30) 


for  the  first  time  step.  By  the  steps  shown  in  Appendix  C,  the  following  is  obtained. 


”]kt  = 


sin(At) 


(N-l)  At 


y%i 

-r^ 

^jkl 

0 


(iV  +  1)  1  —  cos(A<) 
(N-l)  At 


y%i 

0 


The  first  term  on  the  right  side  of  this  equation  is  a  consequence  of  the  difference,  (Wg)y^i 
—  (»p9^2.  appearing  in  Eq.  (30).  Note  that  even  in  the  limit  as  At  —  0,  this  term  does  not 
vanish,  and  the  boundary  velocity  is  given  by 


®}fci 


(N-l) 


y%i 

^jk\ 

0 


Thus,  this  method  creates  nonzero  velocity  components  at  the  wall  which  are  propagated 
into  the  field.  Based  on  the  correspondence  between  terms  in  the  equations  above,  the 
disturbance  is  related  directly  to  the  fact  that  the  grid  velocity  is  evaluated  at  more  than  one 
point  in  Eq.  (25).  On  the  other  hand,  the  magnitude  of  the  initial  disturbance  is  diminished 
with  grid  refinement,  i.e.,  as  AT  —  oo.  Nevertheless,  since  Method  2  generates  this  disturbance 
in  the  quiescent  flow,  it  is  judged  to  be  unacceptable. 

4,5  ANALYSIS  OF  METHOD  3 

Now  consider  the  application  of  Method  3  to  Problem  1  and  then  Problem  2.  For  Problem 
1,  the  curvilinear  coordinate  system  is  time-independent,  so  the  grid  velocity  is  zero.  Also, 
since  Vf  =  e^,  Eq,  (26)  becomes 
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u 

1 

u 

0 

■  0  ‘ 

V 

V 

- 

0 

w 

jkl 

w 

jk2 

1 

0  0  1 


U 

0 

u 

V 

— 

0 

w 

jk2 

0 

for  the  first  time  step.  Note  that  the  last  equality  follows  since  it  is  assumed  that  the  flow 
is  two-dimensional  with  the  v  component  zero.  Since  Method  3  does  not  spuriously  generate 
a  V  component  in  the  field,  it  is  judged  acceptable  in  terms  of  Problem  1. 

Now  consider  the  application  of  Method  3  to  Problem  2.  Since  the  fluid  velocity  is  assumed 
to  be  zero  everywhere  initially,  the  method  becomes 


^  (31) 

for  the  first  time  step.  By  the  same  steps  as  shown  for  Method  2,  the  following  is  obtained. 


,  1  —  cos(At) 

= 


At 


y%i 

0 


As  At  —  0,  the  velocity  introduced  at  the  wall  vanishes.  In  fact,  in  the  limit,  the  grid  velocity 
is  computed  exactly  and  =  0.  On  the  basis  of  these  examples.  Method  3  is  judged 
superior  to  Methods  1  and  2. 

5.0  OTHER  SOLID-WALL  CONDITIONS 


After  Method  3  was  derived,  it  was  later  found  to  be  consistent  with  a  scheme  used  m 
Ref.  4.  A  node-centered  version  of  the  method  for  setting  the  boundary  velocity  is  given  by. 


[vCii-oJ^  +  CO)” 


jkl 


(32) 


That  the  condition  on  velocity  is  equivalent  to  that  shown  in  Eq.  (26)  for  Method  3  can  be 
seen  as  follows.  By  Eqs.  (12),  (2),  and  (5), 

^r) 

Vx  %  Vz  VC  Vv  y< 

Cx  Cy  Cz 


-vc-w,  =  vc^[j-ifJ  =  [o  0  i] 
=  e, .  [jj-^\  €<  =  c,  •  c,  =  Cf 


1 - 

Ct  _ 
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Thus,  the  right  side  of  Eq.  (26)  is 

{vc’h  ■ .?«  -  vcji,  • 

=  -  ]|^y  {vcji  ■  +  (C.6i} 

which  is  on  the  right  side  of  Eq.  (32). 

Next,  the  following  solid-wall  boundary  condition  for  pressure  is  imposed  in  Ref.  5, 

dp 

^  •  »• 

Here,  h  is  an  outwardly  directed  unit  normal  vector  at  the  body  surface,  and  d/dn  denotes 
a  directional  derivative  in  the  direction  of  n.  Also,  represents  the  grid  acceleration.  The 
basis  of  this  boundary  condition  can  be  understood  as  follows.  First,  the  conservation  of 
momentum  can  be  written  as 


dv 


— Vp. 


Here,  dv/dt  denotes  the  total  derivative  of  fluid  velocity  which  is,  of  course,  equal  to  the 
fluid  acceleration,  a.  Now,  in  addition  to  imposing  Eq.  (11)  so  that  the  normal  component 
of  fluid  velocity  and  grid  velocity  are  identical,  one  may  naturally  require  the  same  condition 
for  the  fluid  and  grid  accelerations, 

»  •  («  —  Ug)  =  0. 


Finally,  when  this  equation  is  combined  with  the  normal  component  of  the  preceding  one, 
the  condition  in  Eq.  (33)  is  obtained. 

6.0  CONCLUDING  REMARKS 


This  work  emphasizes  that  subtle  and  erroneous  effects  can  result  from  the  use  of 
apparently  natural  numerical  methods.  It  also  underscores  the  importance  of  developing  simple 
tests  for  these  computational  methods  to  elucidate  the  possible  consequences  of  their  use. 
In  this  report,  simple  model  problems  were  used  to  test  three  methods  for  imposing  discrete 
solid-wall  boundary  conditions  for  the  numerical  computation  of  inviscid  flows.  The  methods 
considered  were  applied  in  stationary  and  moving  grid  systems.  It  was  shown  that  two  of 
these  methods  are  only  provisionally  correct.  In  particular,  a  well-known  method  involving 
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contravariant  velocity  components  was  shown  to  provide  correct  solid-wall  boundary  values 
only  if  grid  lines  emanating  from  the  solid  boundary  are  normally  oriented  at  the  wall. 
Specifically,  this  report  shows  that  without  the  orthogonality  constraint,  the  method  corrupts 
the  tangential  component  of  the  velocity  extrapolated  from  the  field.  Also,  another  more 
direct  method  was  shown  to  provide  correct  solid-wall  boundary  values  only  if  the  grid  velocity 
is  uniform,  as  is  the  case  when  the  frame  of  the  moving  body  is  inertial.  Otherwise,  the 
differences  between  grid  velocities,  evaluated  at  different  points  for  the  method,  disturb  the 
boundary  values.  Finally,  a  correct  procedure  was  posed  to  avoid  the  problems  of  the  other 
two  methods,  and  this  procedure  was  shown  to  pass  tests  developed  herein. 
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APPENDIX  A 

DERIVATION  OF  ORTHOGONAL  TRANSFORMATION 

In  this  appendix  it  is  shown  that  if  a  system  is  in  rigid  motion,  i.e., 

Ilst(t)  —  ®a(t)||  =  ||S(,(0)  —  *a(0)||  (A-1) 

where  sjf)  and  represent  Cartesian  components,  relative  to  the  stationary  frame,  of 
any  two  points  that  are  fixed  in  the  moving  frame,  then  there  is  an  orthogonal  matrix,  L(t),  i.e., 

L-\t)  =  L^{t)  (A-2) 


such  that 


[S6(t)  -  So(t)]  =  L{t)  [S6(0)  -  »a(0)] .  (A-3) 

For  this,  define  a  vector-value  function,!,  by 

£(®(0)  -  ®o(0),t)  =  s{t)  -  so{t).  (A-4) 

where  x{f)  and  ^^it)  represent  Cartesian  components,  relative  to  the  stationary  frame,  of  two 
points  that  are  fixed  in  the  moving  frame.  In  other  words,  I  performs  a  rigid  motion  on  the 
vector,  s(0)  -®q(0),  and  produces  x(t)  -  after  time  t.  To  achieve  the  result  of  Eq.  (A-3), 
fwill  be  shown  to  be  linear  in  its  first  argument  and  representable  in  terms  of  an  orthogonal 
matrix. 

First,  it  is  shown  that!  preserves  scalar  products,  i.e.,  for  any  vectors,  x  and  y. 


i{x,t)- l{y,t)  =  X -y.  (A-5) 

For  this,  x{t)  and  y{i)  represent  points  fixed  in  the  moving  frame  such  that 

*(0)  —  xo(0)  =  X  and  f  (0)  —  xo(0)  =  f.  (A-6) 

Then,  using  the  identity 

2®  •  w  =  ||®||^  +  ll®ll^  ~  II®  ~  ®||^ 


25 


AEDC-TR-93-23 


and  the  definition  of  t  gives 

€(®(0)  -  *o(0),t)  •€(y(0)  -  *o(0),t)  =  [®(t)  -  xo(t)]  •  [y(t)  -  *o(0] 
=  i||*(t)  -  *o(t)||H  i||,(t)  -  xoWlP  -  i||*(t)  - 


By  Eq.  (A-1), 

«(x(0)  -  xo(0),t)  •  «(s(0)  -  xo(0),t)  =  i||x(0)  -  xo(0)||2  +  i||9(0)  -  xo(0)||2 
-^ll*(0)  -  !r(0)|p  =  [x(0)  -  xo(0)] .  [,(0)  -  xo(0)] . 

Equation  (A-5)  now  follows  from  the  above  equations. 

This  fact  is  now  used  to  show  that 

\\e{a3;  +  I3y,t)  -  at{x,t)  -  I3e{y,t)\f  =  0 

and  hence  that  €  is  linear  in  its  first  argument.  Using  Eq.  (A-5), 

||£(q:x  +  /3y,t)  -  al{x,t)  -  /?€(y,t)|p 
=  i{ax  +  /3y,t)  •  £(ax  -|-  py,t)  +  a‘i{x,t)  •i{x,t)  +  l3‘£{y,t)  •£(§,<) 

-2ai(ax  -1-  0y,t)  •  i{x,t)  —  2fi£{ax  +  ^§,t)  ■  £{y,t)  —  2a(}£{x,t)  •  £(y,t) 

=  {ax  +  I3y)  ■  {ax  -f  /3y)  -F  a^x  ■  x  +  0^y  ■  y 

—2a  {ax  -f-  Py)  •  s  —  2/?  {ax  -f  py)  ■  y  —  2a/?®  •  y 
=  \\{ax  +  Py)  -  ax  -  Py\f  =  0. 

Thus,  I  is  linear  in  its  first  argument  and  can  be  written  in  terms  of  a  matrix,  L{t),  as 

i{x,t)  =  L{t)x.  (A-7) 
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To  see  that  Lit)  must  be  orthogonal,  note  that  by  Eqs.  (A-5)  and  (A-7) 

X  •  y  =  £{x,t)  •  £iy,t)  =  x^X^(t)i/(t)§. 

Since  *  and  y  are  arbitrary,  they  can  be  chosen  as  i  =  e.  and  f  =  e^.  to  obtain 

ejL'^it)Lit)eJ  = 

In  other  words,  Eq.  (A-2)  holds.  Finally, 

Lit)  [S6(0)  -  Sa(0)]  =  Lit)  {[ar(,(0)  -  xo(0)]  -  [sa(0)  -  3Co(0)]} 

=  ^(xt(0)  -  Xo(0))0  —  ^(®a(0)  -  Xo(0),t) 

=  [Xbit)  -  Xo(t)]  -  [Xa(t)  -  Xo(t)] 

—  * 

Thus,  Eq.  (A-3)  is  obtained. 
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APPENDIX  B 

DERIVATION  OF  CONTINUUM  CONDITION 

In  this  appendix,  the  continuum  condition  of  Eq.  (10), 

n-{v-  Wrn)  =  0  (B-1) 

is  derived  from  the  more  primitive  condition, 

«  •  «  =  0  (B-2) 

which  applies  in  the  moving  frame.  Specifically,  n  represents  the  Cartesian  components,  relative 
to  the  moving  frame,  of  an  outwardly  directed  unit  normal  vector  at  the  body  surface.  Also, 
V  represents  the  Cartesian  components  of  fluid  velocity  in  the  moving  frame. 

To  establish  the  relationship  between  the  components  of  Eq.  (B-1)  in  the  stationary  frame 
and  the  components  of  Eq.  (B-2)  in  the  moving  frame,  the  two  coordinate  systems  are  related 
as  follows.  Specifically,  let  *  and  i  represent  the  Cartesian  components,  relative  to  the 
stationary  and  moving  frames,  respectively,  of  any  point  at  time  t.  Also,  suppose 
represents  the  Cartesian  components,  relative  to  the  stationary  frame,  of  a  reference  point 
used  as  the  origin  in  the  moving  frame.  Assume  that  the  two  coordinate  systems  coincide 
at  t  =  0  and  thus  ^^^(0)  =  0.  Then  by  Eq.  (9),  L\t){x  represents  an  initial  vector 

which  from  time,  f  =  0,  evolves  to  the  vector  L{t)L\t)[x  =  x  -x^it)  after  time  t. 

Furthermore,  the  vector  remains  in  rigid  motion  with  respect  to  the  moving  frame  and  hence 
has  unchanging  Cartesian  components,®,  relative  to  the  moving  frame.  Therefore,  since  the 
two  coordinate  systems  coincide  initially, 

X  =  L"{t)  [x  -  ®o(0]  •  (B-3) 


Also,  by  Eq.  (8), 

X  —  xo(t)  =  L{t)x.  (B-4) 

Equations  (B-3)  and  (B-4)  define  the  relation  between  the  stationary  and  moving  frame 
Cartesian  systems. 

Now,  the  relation  between  n  of  Eq.  (6)  and  »  will  be  determined.  For  this,  let  xjf) 
represent  the  Cartesian  components,  relative  to  the  stationary  frame,  of  a  point  on  the 
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boundary.  Then  define  =  xjj)  +  h{t)  to  represent  a  point  off  the  wall  at  the  end  of 
the  unit  normal  vector,  n.  Since  at  t  =  0,  ®^(0)  —  *^(0)  =  n(0)  =  n,  Eq.  (9)  gives 

h{t)  =  [xb(t)  -  *a(0]  =  -^(0  [®6(0)  -  Xa(0)]  =  L{t)h{0)  =  L{t)n.  (B-5) 

Also,  by  Eq.  (8), 

n  =  L'^{t)n(t).  (B-6) 


Next,  V  will  be  related  to  ®  and  Specifically,  it  will  be  shown  that 

w(t)  =  [®(^)  ®m(0]  •  (B"7) 


For  this,  let  Sj(t)  and  ij(t)  represent  the  Cartesian  components,  relative  to  the  stationary  and 
moving  frames,  respectively,  of  a  fluid  particle  trajectory.  Also,  let  and  represent 
the  Cartesian  components,  relative  to  the  stationary  and  moving  frames,  respectively,  of  the 
trajectory  of  a  material  particle  at  rest  in  the  moving  frame.  Then  define 


V  = 


®m 


dXm 

~dr' 


fi  ^ 

“  ~dr  ~ 


dxQ 


Equation  (B-7)  will  now  be  established  for  an  arbitrary  point  in  space-time  where  the  fluid 
particle  and  material  particle  trajectories  intersect.  Specifically,  it  is  assumed  that  there  is 
a  time,  t*,  when 

-  Zjn{t*)  and  Xf{t*)  =  Xm-  (B-8) 


Now,  by  Eq.  (B-3), 


Xf{t)  =  L'^it)[xf{t)-xoit)\. 


Differentiating  this  gives 


jtT 

vit)  =  -  so(i)]  +  L^{t)  WO  -  »o(0] 


By  Eq.  (B-4), 


Xf{t)  -  ®o(0  = 
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Combining  the  last  two  equations  gives 


firT 

v{t)  =  —{t)L{t)xf{t)  +  L^{t)  [v{t)  -  Oo(<)] 


(B-9) 


To  bridge  the  gap  between  Eq.  (B-9)  and  Eq.  (B-7),  the  material  particle  velocity  will  be 
expressed  in  terms  of  the  items  appearing  on  the  right  side  of  Eq.  (B-9).  By  Eq.  (B-4), 


By  Eq.  (B-8), 


*7n(0  *o(0  —  L{i)Xm{fy  —  L(t)Xm. 

~  ®o(0  ~  )• 


After  differentiation  this  becomes 


dL. 


®m(0  -  »o(0  =  —it)xj{t*). 


Multiplying  this  by  L\t)  and  evaluating  the  result  at  t  =  t*  gives 


'J'  y  dL  , 


{t*)  K(r )  -  voin]  = 


To  relate  the  right  side  here  with  the  term  involving  in  Eq.  (B-9),  the  following  relation 
is  used, 


dL'^ 


dL. 


it  = 


Since  L{t)  is  orthogonal,  this  follows  by  differentiating  L^{t)L{t)  =  I.  From  the  last  two 
equations  it  follows  that 


-  i^(t*)  [wm(<*)  -  vo{t*)] 


dlT 

di 


(B-10) 
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Evaluating  Eq.  (B-9)  at  /  =  /*  and  combining  the  result  with  Eq.  (B-10)  gives 

«(r)  =  -i^(r)K(r)-i,o(t*)]  +  x^(t*)[»(r)-oo(f*)] 

Since  the  intersection  point  of  Eq.  (B-8)  is  arbitrary,  Eq.  (B-7)  follows. 

Now,  combining  Eqs.  (B-2),  (B-6),  and  (B-7)  gives  the  continuum  condition  on  the  fluid 
velocity  expressed  in  terms  of  Cartesian  components  in  the  stationary  frame, 

=  rJ LL^  [o  -  Vm]  =  n  •  [»  -  Wm] 
which  is  the  condition  of  Eq.  (B-1). 
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APPENDIX  C 

DETAILED  CALCULATIONS  FOR  ANALYSIS  OF  METHOD  2 

The  purpose  of  this  appendix  is  to  provide  detailed  calculations  to  demonstrate  the  result 
of  applying  Method  2  in  Eq.  (25)  to  Problem  2  defined  in  Subsection  4.2.  Since  the  fluid 
velocity  is  assumed  to  be  zero  everywhere  initially,  the  method  becomes 


(C-l) 


for  the  first  time  step.  To  compute  the  grid  velocity  for  this  by  Eq.  (24),  the  orthogonal 
transformation,  L{t),  must  be  computed.  For  this,  take  the  reference  point,  x^it)  =  0,  so 
that  =  L(t)xji^f()).  Then  writing  in  terms  of  s^^/O)  defines  Lit)  as  follows.  By 
Eq.  (29), 


Zjklit) 


l  +  iN-2) 
(N-l) 

l  +  iN-2) 
(N-l) 


COS 


t- 


'  cos(t)  cos 


(fc-l)1 
(N  -  1)J 

(fc-1) 


(iV-l)J 


sin(t)  sin 


(iV-l)Jj 


=  cos(t).Tjfc/(0)  -  sin(t)yjfc/(0) 


Vjklii)  = 


l  +  iN-2) 
(N-l) 

l  +  iN-2) 
(N 


sin 


t- 


(fc-1) 


|sin(t) 


(N-l) 
ik 


cos 


1) 


L  (i'z-i) 


+  COS 


sin(t)a;jfc/(0)  +  cos(<)i/jfc/(0) 

U  - 1) 

(N-l) 


Therefore,  defining 


L{t)  = 


cos(t)  —  sin(t)  0 
sin(t)  cos(t)  0 

0  0  1 
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satisfies  the  equation,  Thus,  by  Eq.  (24)  the  grid  velocity  is  given  by 


I- 


cos(t”  —  At)  —  sin(t"  —  At)  0 
sin(t”  —  At)  cos(t"  —  At)  0 
0  0  1 


cos(t”)  sin(t")  0 
—  sin(t")  cos(t”)  0 
0  0  1 


^jkl 


which  simplifies  to 


Now,  since 


1  —  cos(  At)  —  sin(  At)  0 

sin(At)  1  —  cos(At)  0 
0  0  0 


(C-2) 


«o  _ 


COS 


y%i 


^jkl 


=  Sin 


(j  - 1) 


- 

^jk2  - 


(i  - 1) 


(IV -1) 

the  first  two  terms  in  Eq.  (C-1)  are  given  by 

1 


(iV-1) 


At 


1 

At 


1  —  co.s(At)  —  sin(At)  0 
sin(At)  1  —  cos(At)  0 
0  0  0 


(A;-!)- 

0  ^ 

(A;-!)- 

N 

(N-l). 

^Jk2  -  (jV  _  1) 

(N-1). 

(N-1) 

(^-1)1 

0  _  iV  . 

■  (A:-l)l 

N 

(iV-l)J 

yjk2 

-  (^-1). 

(N-l) 

-r*?. 

■  ^jkl 


(C-3) 


{N  -  l)At 


-  ^%2)  (1  -  cos(At))  -  sin(At) 

-  ^'jk2)  sin(At)  +  -  j/9^.2)  (1  -  cos(At)) 

0 


-  cos(At))  -  sm(At) 
sin(At)  +  -  cos(At)) 

0 


(C-4) 


Next,  by  Eq.  (28),  in  Eq.  (C-1)  is  given  by 
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vc^i  = 


(^-1) 


+  {y%i) 


r  r®  1 

^jkl 

r®  1 
^jkl 

J 

y%i 

=  (iV-l) 

y%i 

0 

0 

(C-5) 


since  ^  the  body  surface.  Now,  combining  Eqs.  (C-2)  and  (C-5)  gives 

vcj^i  • 

N  -1 
At 

N  -  1 


’  -  cos(At))  -  y9,^^sm{At)  ' 

,  y%i  0  ] 

x]k2  sin(At)  +  -  cos(At)) 

^}k2 

At 


{sin(At)  [x%2yjki  -  ^%iy%2]  +  (1  -  cos(At))  [x%^x%2  +  y%iyjk2] }  ■ 


Using  Eq.  (C-3)  to  simplify  this  gives 


xo  _^rl-cos(At) 
^Cjkl  •  i^9)jk2  -  ^ 


Thus,  according  to  this  equation  and  Eq.  (C-5),  the  last  term  in  Eq.  (C-1)  is. 


llVC^alP 


'^Cjkl  '  {^g)jk2  ~ 


N  1  —  cos(At) 
{N  -  1)  At 


®ifci 

ym 

0 


Finally,  combining  Eqs.  (C-4)  and  (C-6)  for  Eq.  (C-1)  gives 


1  1  sin(At) 

1 

i 

1 _ 

(iV  -1- 1)  1  —  cos(A/) 

®ifci  1 
.0 

At 

-^jk\ 

0 

(tV-1)  At 

- 1 

tH 

/r\ 


This  provides  the  required  result  for  Subsection  4.4. 
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NOMENCLATURE 

a  .  Fluid  acceleration 

.  Grid  acceleration 

V .  Space-time  gradient  operator  defined  in  Eq.  (14) 

*>>>  ®z»  •  Cartesian  basis  for  space-time 

j,  k,  I  .  Grid  point  indices  related  to  curvilinear  coordinates  according  to  Eq.  (1) 

J .  Jacobian  of  the  transformation  from  *  to  C  defined  in  Eq.  (2) 

I .  Vector-valued  function  defined  in  Eq.  (A-4)  for  the  derivation  of  the 

orthogonal  transformation. 


L(?)  .  Orthogonal  transformation  defined  in  Eq.  (9) 

M  .  Constant  matrix  defined  in  Eq.  (19) 

n  .  Time  level  index  related  to  curvilinear  coordinates  according  to  Eq.  (1) 

N .  Grid  refinement  parameter  defined  following  Eq.  (28) 

n,  »  .  Cartesian  components,  relative  to  the  stationary  and  moving  frames, 

respectively,  of  a  vector  normal  to  the  solid-body  surface 

p  .  Pressure 

t  .  Time 

t* .  Time  at  which  the  fluid  particle  and  material  particle  trajectories  intersect 

as  shown  in  Eq.  (B-8) 

M,  V,  w .  Cartesian  components  of  fluid  velocity  relative  to  the  stationary  frame 

U,  V,  W,  S  .  Contravariant  components  of  V 
li .  Generic  space-time  vector  used  in  Subsection  3.1 


35 


AEDC-TR-93-23 


V .  Temporal  derivative  of  IV 

h .  Cartesian  components  of  fluid  velocity,  relative  to  the  stationary  and  moving 

frames,  respectively 

Wq  .  Cartesian  components,  relative  to  the  stationary  frame,  of  the  velocity  of 

a  reference  point  whose  trajectory  is  given  by  in  the  stationary  frame 

.  Cartesian  components,  relative  to  the  stationary  frame,  of  the  velocity  of 

a  point  with  fixed  curvilinear  coordinates,  and  trajectory,  in  the 
stationary  frame 

.  Cartesian  components,  relative  to  the  stationary  and  moving  frames, 

respectively,  of  the  velocity  of  a  material  particle  whose  trajectory  has 
Cartesian  components,  and  x^,  in  the  stationary  and  moving  frames, 
respectively 


W  .  World  line  of  a  fluid  particle,  used  in  Subsection  3.1 

X,  y,  z  .  Stationary  frame  Cartesian  coordinates 


Xj,  Xj,  Xj  ...  Same  as  x,  y,  z,  respectively. 

I,  i .  Cartesian  components  of  a  point  relative  to  the  stationary  and  moving 

frames,  respectively.  Also,  x  represents  [x,  y,  z]^ 

Xq  .  Cartesian  components,  relative  to  the  stationary  frame,  of  the  trajectory 

of  a  reference  point  with  fixed  curvilinear  coordinates 

®y.,  .  Cartesian  components,  relative  to  the  stationary  and  moving  frames, 

respectively,  of  the  trajectory  of  a  fluid  particle 

X  .  Cartesian  components,  relative  to  the  stationary  frame,  of  the  trajectory 

O 

of  a  point  with  fixed  curvilinear  coordinates, 

x^,  .  Cartesian  components,  relative  to  the  stationary  and  moving  frames, 

respectively,  of  the  trajectory  of  a  material  particle 

X,  f .  Arbitrary  vectors  used  in  Eqs.  (A-5)  and  (A-6) 
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X . 

T  . 

^0-  . 

Af  . 

r;,  f 

^2’  ^3 

^  . 

. 


T  .  .  . 


(.f 


(•).. 


(•)... 


Space-time  position  vector 

Parameter  used  in  the  definition  of  the  curvilinear  coordinate  system  of 
Problem  1  as  shown  in  Eq.  (27) 

Bounds  on  det{J}  as  shown  in  Eq.  (3) 

Kronecker  delta,  where  =  0  if  i  ^  j,  and  b-j  =  1  if  /  =  y 

Time  step 

Curvilinear  coordinates  established  around  the  moving  body. 

Same  as  t/,  f. 

Represents  [|,  j;,  z]^ 

Fixed  curvilinear  coordinates  of  a  point  with  trajectory  having  Cartesian 
components,  s^,  relative  to  the  stationary  frame 

Fixed  curvilinear  coordinates  of  a  point  with  trajectory  having  Cartesian 
components,  relative  to  the  stationary  frame 

Curvilinear  coordinate  to  complement  j;,  f  in  space-time. 

Subscripts,  j,  k,  /,  and  superscript,  n,  indicate  space-time  grid  location 

Superscript,  T,  indicates  transposition  of  a  matrix 

Subscripts,  indicate  differentiation  with  respect  to  the  respective 

variables 
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